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The linear instability of a beam tensioned by its own weight is considered. It is shown that for long 
beams, in the sense of an adequate dimensionless parameter, the characteristics of the instability 
caused by a follower force do not depend on the length. The asymptotic regime significantly differs 
from that of short beams: flutter prevails for all types of follower loads, and flutter is localized at 
the edge of the beam. An approximate solution using matched assymptotic expansion is proposed 
for the case of a semi-infinite beam. Using a local criterion based on the stability of waves, the 
characteristics of this regime as well as its range of application can be well predicted. These results 
are finally discussed in relation with cases of flow-induced instabilities of slender structures. 
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I. INTRODUCTION 

The linear stability of a beam under the action of a follower force exerted at one of its extremities has been the 
subject of intensive research, as can be seen from the extensive review in m- The interest in this problem lies in 
its potential applications and also in the large variety of fundamental topics of mechanics involved in the solution 
method. Practical examples of direct application in the field of fluid-structure interactions are numerous: fluid- 
conveying pipes, plates subjected to axial flow or towed cylindrical bodies are modelled by equations that are similar, 
though not identical, to those of a beam under a compressive follower force [I113- More generally, follower forces 
have been extensivelly discussed in the literature, including in terms of their physical reality, see Q for a full review. 
The case of a cantilevered beam of finite length under a partial follower force is well documented, see 0 , with many 
results on the effect of characteristics of the beam or of the load on the critical load that causes instability, and on 
the nature of the instability, be it divergence (buckling) or flutter. 

We seek here to establish the characteristics of instability of a beam in the case where its length is much larger 
than the region where an unstable motion will develop. This arises when a constant load, such as gravity acting on 
a vertically hanging beam, produces a tension that increases along the beam, from zero at the lower free end to a 
maximum at the upper fixed end. The increasing tension induces a corresponding increasing stiffness. Motion is then 
confined to the lower end, corresponding to edge flutter. This has been studied both experimentally and numerically 
in three of the problems of flow-induced vibrations mentioned above: hanging fluid-conveying pipes 0,li3> hanging 
ribbons under axial flow [T^ and towed cylinders under axial flow S0. In all these systems it was observed that 
there exists a limit state in which the length does not affect the stability. This limit state is found for length larger 
than a limit value given by simple considerations on the local stability of bending waves SB- 

The objective of this paper is to establish similar results on the generic case of a beam under a partially follower 
force (or subtangential force), tensioned by a load such as gravity. 

In Section |TT] we shall recall the equations of motion and the possible choices of dimensionless variables. The effect 
of the beam length on stability is analysed in Section Hill using numerical computations of the eigenmodes. In Section 
we address the particular case of a semi-infinite beam, using various types of modelling. The application of the results 
given in the paper to problems of flow-induced instabilities is discussed in Section IVl 


II. EQUATIONS OF MOTION 

We consider a vertical beam of length L loaded by its own weight, FiglT] 
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FIG. 1: Semi-infinite hanging beam with a partially follower force. 


A partial follower force is applied at its lower end m- The linear equation governing the in-plane lateral deflection 
Y{Z,T) reads 






dY' 

{P-mg Z) — 


d^Y 


( 1 ) 


where El is the flexural rigidity, P is the load, g is gravity and m is the mass per unit length. No daiming is 
considered here, though it is known to significantly influence some aspects of the problem, see for instance in Q. The 
boundary conditions at the lower end, Z — 0 are 


q2y q^y dY 

= 0 , + (1 - V)P^{0) = 0 , 


( 2 ) 


where 77 is a coefficient that expresses the part of the loading that follows the beam slope; hence 77 = 0 corresponds to 
a non-follower force and 77 = 1 to a pure follower force. At the upper end, we assume a clamped condition. 


. ^ dY, , 

Using the length of the beam, L, as a reference we define the following dimensionless variables: 
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Then Eq.([T|) may be re-written in dimensionless form as 
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with the boundary conditions 

d^y 


7^,2(0)=0; §(0) + (i-^>|(0) = 0; 2/(1) =0; |(1) = 0. 


(3) 


(4) 


(5) 


( 6 ) 


This dimensionless set of equations is adequate to analyse the effect of 7 on the critical loading, but only if the length 
L is kept constant. As we need to vary the length L, it is necessary to define a new set of dimensionless parameters. 
We now use a length scale defined by the ratio of the two stiffnesses of the beam, namely the flexural rigidity and the 
stiffness related to tension induced by gravity Q, 
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which will be referred to as the gravity length. Note that this length scales like the critical length that causes buckling 
of the beam under its own weight. Using Lg as the reference length to define dimensionless variables we have 


X 

where Lg has been substituted 

x = zi, T = U\q=^, l = 

Eq. (U) now becomes 
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for L. These are related to the previous set of variables by 


dx^ 


dx 


{q-x) 


dx 



( 8 ) 

(9) 


( 10 ) 


and the corresponding boundary conditions are 




(11) 


III. STABILITY OF FINITE BEAMS 
A. Solution for short beams 

For short beams, in the sense where <C 1, the critical load qc may be derived by considering that 7 <C 1 so that 
Eq. (0 becomes 


dz dP 


( 12 ) 


The corresponding stability diagram of Pc versus rj is that of the generalized Beck’s column m and is shown in 
FigUJa). For rj < 0.5 instability arises in the form of a divergence, whereas flutter prevails for rj > 0.5. The 
corresponding frequency at the flutter limit is shown in Figl2l(b). The critical load and frequency at flutter expressed 
in the dimensionless variables of Eq. m are 


QcivJ) 


Pcjy) 

e ’ 


Wc(?7,-^) = 


Uc(??) 


(13) 


B. Effect of length 

We now investigate the effect of the length I on the stability threshold and the nature of the instability. The 
numerical method utilized to this end is the finite difference method, as in [a. It allows us to compute the eigenvalues 
uj and the eigenvectors such that 


+ [(9 “ x )^]' ~ = 0 


(14) 


with 


(p"(0) = (/?(^^(0) + (1 - v)qp'{0) = ip{£) = ip'{£) = 0. (15) 

The critical load and the nature of the corresponding instability can be determined from the evolution of the real and 
imaginary parts of the frequency. The beam is unstable if the frequency w has a negative imaginary part. If the real 
part is zero, the instability results in an exponential growth in time of the deformation , without oscillation and the 
instability is of the divergence type. If the real part is non zero, the instability results in an exponential growth of 
oscillations, and the instability is of the flutter type. 

FigElshow the evolution of the critical load with length, for several values of rj of particular interest. 







4 




FIG. 2: Critical load, (a), and frequency at instability, (b), for short beams. Dashed vertical lines are guides for the eyes that 
do not correspond to calculations. 


For 77 = 0, Figl3{a), which is the case of a non-follower force, the critical load decreases steeply with length, up to 
about £ = 5 where it reaches a limit value and does not change when the length is increased further. Instability is for 
all lengths of the diver genc e type. On the same figure are shown computations from [I^ and [I^, for intermediate 
values of the length. In , the static instability of a towed beam has been considered, which yields equations similar 
to that used here. Details of the equivalence are discussed in the last section of this paper, but suffice to say here 
that we may use the results of Fig. 4 of 0 with the change of variables I = and q — . It is seen that their 

results and those of this paper are in very good agreement. For 77 = 0.2, FiglSjb), the critical load for divergence 
also decreases and reaches a limit value, but the flutter threshold becomes lower than that for divergence, for beams 
longer than £ = 5. For 77 = 0.3, Figl^Jc), no divergence is found for £ > 3. Instability is that range of £ is of the flutter 
type. For 77 = 1, FiglSjd), the system loses stability by flutter, similarly to Beck’s column; the critical load for flutter 
decreases with £ until it reaches a limit value. 

From these four cases it can be stated that there exists a limit configuration for all long beams, say £ > 10. For 
these lengths, the type of instability, flutter or divergence, may differ from that observed for short beams at the same 
value of 77 . 

This is further confirmed in the evolutions of the stability diagram qcij\) shown in FigS) At £ = 1 the diagram is 
similar, though not identical, to that for short beams; since £ = 1 we have here qc = Vc and Figs. [2]and[3]^a) can be 
directly compared. At £ = 100, Fig[4l[b) a limit state is almost reached, where stability is always lost by flutter, for 
all values of 77 . 

FiglHshows the mode shape at the critical load, which may be a divergence or a flutter instability depending on the 
values of 77 and £, as shown above. Note that at the instability threshold kp is real, even for flutter instability. In fact, 
in Eq. (HI, when bj is real so is p>. For 77 = 0, Fig[5l^a), the mode shape, here for divergence, becomes independent 
of the length when £ > 3. This is consistent with FiglD^a) where the critical load was found to be stationnary in this 
range. For 77 = 0.3 and 77 = 1, Figs. EJb) and (c), the mode shapes also converge to a constant shape as the length 
is increased. There, the displacement is confined to the lower part of the beam. Note that the instability is of the 
flutter type, except for 77 = 0.3 at £ = I. Yet, the limit mode shape closely resembles that of the divergence instability 
at 77 = 0 . 


IV. STABILITY OF A SEMI-INFINITE BEAM 
A. Boundary conditions 

We now seek to determine the characteristics of the instability when the beam is assumed to be of infinite length in 
the X axis. This is expected to be the solution to which the results of the preceding section converge as £ is increased. 
Whereas the boundary condition at the lower end is unchanged, the condition that the beam is clamped at X = L 
needs to be replaced by the conditions that the displacement remains finite as X goes to infinity and that in this limit 
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(a) 


(b) 



(c) (d) 

FIG. 3: Effect of length on the critical load for divergence and for flutter instabilities, (a) rj = 0,(b) 77 = 0.2, (c) rj = 0.3, 
(d) 77 = 1.. (—), present results for divergence; (- -), present results for flutter; (o), computations by [l^ for 77 = 0; (A), 
computations by [l^ . 


propagating waves only radiate in the direction of increasing X. These conditions in dimensionless form read, 


lim \y\ = 0 ; 

x—^oc 


y dydy 

x^oo Ox Ot 


(16) 


B. Divergence instability 


The case of divergence instability may be analysed by neglecting all time derivatives in the equations, so that (11011 
becomes 


A 

dx^ 


A 

dx 


( \^y 


= 0 


(17) 


and the boundary conditions are only 

d'^y 


dx^ 


(0) = 0; 


0(0) + (1 - 77)77^(0) =0; Jim |j/|=0. 


(18) 
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FIG. 4: Effect of the ratio rj on the critical load for divergence, (—), and for flutter, (- -). (a) Short beam, £ = 1, (b) Long 
beam, I = 100 . 



FIG. 5: Mode shape at the critical load, (a) 77 = 0, (b) 77 = 0.3, (c) 77 = 1. (..), £=!,(- -), £ = 3; £ = 5; (—), £ = 10. 


By defining A = and Z = x — q, this set of equations is identical to that solved by for the case of a towed 
cylinder. Using their solution with our variables, the critical load for divergence qciiq) is found to satisfy 


where Ai is the Airy function 


a- 


Ai'(-gc) - qd 


1 

3 


Ai(s)ds 


= 0 , 


(19) 


C. Flutter instability 

We give here an approximate solution for the critical load for flutter of a semi-infinite beam, which is an extension 
from the work of |10l| for a purely follower force, 77 = 1. The approach used in the matched assymptotic expansions 
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follows that proposed by [l9j | for the dynamics of vibrating strings tensionned by gravity. 

We consider Eq. with the boundary conditions m at the bottom end, and dT51) for a semi-infinite beam. We 
seek harmonic solutions of this set of equations, of the form 

y{x,t) = Keal[ip{x)e^‘^*] (20) 

as a function of the follower force magnitude q. Instability will be associated with Im(a;) < 0. The variation in space 
of the solution satisfies 


+ \{q-x)(f']'- uj'^lP = Q, (21) 

with the boundary conditions 

(/j(2)(0) = + (1 _ y)qq,'{Q) = 0 (22) 

where ()' denotes derivation with respect to x. The condition for radiating waves reads 

lim (Real[(/7'(a:)e“‘]Real[ia;(/3(a;)e*‘^‘]) < 0. (23) 

Let us consider first the lower part of the beam where x is of the order of q. In this “inner” domain of length g, we 
may derive an approximate solution ipi (x) by simply lumping all its inertia at the lower end of the beam, so that the 
boundary conditions at x = 0 read now 

(^p)(0) = 0;</jf^(0) -b (1 - ??)Q‘b’U0) = 9wV»(0), (24) 

and Eq. (1^ reduces to 

+ [(9 - x)(p'i\' = 0. (25) 


The corresponding solution reads 


(pi{x) = a + / [5Ai(s — g) -b cBi(s — g) -b dGi(s — g)] ds (26) 

Jo 

where Ai, Bi and Gi are the Airy functions, a, 6, c, d being four coefficients. For the sake of clarity we now use the 
notation F{x) = F(s)ds and F for F{—q) unless otherwise specified, where F could be any of the Airy functions 
above. 

In terms of the coefficients a, b, c, d the boundary conditions at the lower end are 

6Ai' -b cBi' -b dGi' = 0; (27) 

6Ai" -b cBi" -b dGi" -b (1 — r])q{hAi + cBi -b dGi) = aquj^. (28) 

Conversely, in the upper part of the beam, that is where x ^ q, we rescale the equations by using new variables in 
space and time, namely 


X = xe; r= (29) 

Ve 

where e = Lg/A, the length A being an arbitrary large scale. Using this set of variables, Eq. (|2T]) reads 

+ [{e^q — £x)v'\' ~ = Oj (30) 

where {)' denotes now differentiation with respect to y. This allows us to derive the equation at the leading order in 
the “outer” domain. 


{-X’f'e)' - = 0, 


the solution of which, in terms of Bessel functions, is (see 


also in Q) 


(31) 


ipeix) = + l3Yo{2ry/x)- 


(32) 


The radiation condition, Eq. (|23)) . implies that 


a — i/3 = 0. 


(33) 


We may now match the inner and outer solutions by considering their respective limits Equating these two 

limits yields three conditions for the coefficients defining the inner and outer solutions, namely 


c = 0; d = /3; a + b 


- Ai 


+ d 


27 + In 3 
Stt 


-Gi 


= a + /3— (Inw + 7 ) 

TT 


(34) 


where 7 is Euler’s constant. The set of boundary conditions (071) . (1751) . (1551) . with the three matching conditions (IMl) 
are six linear equations between the six coefficients a, /3, a, &, c, d. This defines an implicit relationship between 
the two parameters uj and q, so that there exists a non-trivial solution. It reads 


- T]q[Gi'Ai - GiAi'] 


—Ai' + quj'^ 


/I p-' /ln3 — 47 — 3f7r — 61na; 

Vs ■ V “ V 3)^ 




= 0 . 


(35) 


The particular value of the flutter instability threshold qc may be directly derived by assuming that w is real in 
Eg. lfM]) . As all functions of q in (1351) are real this implies that 

Ai’{-q,) = 0; Q, ~ 1.02. (36) 

The real part of the frequency at the instability threshold is then derived using ([35] ) and ([36| ) as 


2 ^ Ai(-ge) 

^ {1/3) - Ai{-qc) 


(37) 


At 7 / = 0, the solution of [T^ for divergence is recovered, both in terms of critical load and frequency. Using this 
approximate solution, flutter is found to exist for all values of rj, and the corresponding critical load does not vary 
with 7 . In Fig. | 6 Ka) this is compared with the numerical solution for very long beams, £ = 100. The approximate 
solution, equation (1551) . predicts a lower bound, qc = 10.2, equal to that for rj = 0. The flutter frequency is very well 
predicted, as may be seen in FiglSKb). At the instability threshold the modal shape ip is real and reduces to 


r 1 

(p{x) = / Ai{s — p)ds — —. (38) 

Jo 3 

This is compared, FigjTl with the computed mode shapes at £ = 10. For ry = 0, FigjTKa), Eq. (1551) is actually the exact 
solution for a semi-infinite beam, which derives from the results of [l^ . For 77 = 0.3 and even 77 = 1, it is remarkable 
that Eq. (1551) stills gives a very good approximation of the mode shape at the flutter threshold. Figs. mb) and (c). 


D. A model based on local wave stability 

Following 0, we consider now a criterion based on the local characteristics of waves, at a given position x. From 
Eq. (ITUl) . the dispersion relation is 

— {q — x)k^ — = 0 (39) 

where k is the wavenumber and lu is the frequency of the wave. All points such that x < q bear unstable waves. 
Beams of length £ larger than q are expected to have a behaviour not affected by length as waves are damped above 
X = q 

In dimensional variables this allows us to define a the length Ly for neutral stability, 

Ln = -, (40) 

mg 

above which the medium only bears stable waves. In Eq. (1) this is the point where the local tension goes through 
zero. Gonsidering that the part of the beam above this, Z > Ljq, plays a negligible role in the instability, we may 
approximate the critical load for the semi-infinite beam by that for a finite beam of length Ljv: a beam of length Ljv, 
without gravity, has a dimensional critical load 
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FIG. 6: Critical load q^, (a), and frequency Wc of flutter,(b), of a long beam. (—), numerical results for I = 100; (-.-) solution 
for the semi-infinite beam using matched assymptotic expansions, Eqs. (1361) and (Eij. 



FIG. 7: Mode shape at the critical load for a long beam, (a) rj = 0, (b) — 0.3, (c) r; = 1. (—), numerical results for i = 10; 

(o), solution for the semi-infinite beam using matched assymptotic expansions, Eq. HSU- 


where Pc is the dimensionless critical load for short beams. As the length Ln varies with the load P, this results in 


Pc = EI [mg f Pe¬ 


rn 


Using now dimensionless variables pertaining to the case of long beams, this reads 


qc=P 


1 / 3 . 

c ’ 


(43) 


similarly, we have 


2/3 ■ 
Pc 


W, 


(44) 
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FIG. 8: Models for the critical load based on local wave stability, (a) critical load (—) numerical results for I = 100; (- -), 
model based on local wave stability, Eq. (1431) and (1441) . (b) Effect of length r) = 1. (—), computations of Section [III Bl (- -), 
short beam model; (...), transition line; long beam model. 


In Fig. [D^a) this is plotted for comparison with numerical results for the very long beam. The order of magnitude 
of the critical load is well recovered. Yet, divergence is predicted for 77 < 0.5, instead of flutter. Note also that the 
critical load at 77 = 1 is very well predicted. Figure Hl^b) summarizes this approach: the limit line £ = g is found to 
be a good approximation of the transition between the behaviour of short beams, that depend on length, and that 
of long beams, that does not. For short beams the approximation of Eq. m applies. For long beams Eq. (l4^ is a 
good approximation. 


V. DISCUSSION 
A. Domains of instability 

In the results presented in this paper two regions can be indentified in a beam tensionned by a load such as gravity. 
In the part near the free end, a; 1 in our dimensionless variables, that \s Z Lg = {El/MgY^^, the stiffness that 
opposes lateral displacement is dominantly that of the flexural rigidity. Conversely for x >> 1 the stiffness is due to 
the tension that results from gravity. 

For a given beam, if its length L is much smaller than Lg, or <C 1, gravity effects can be neglected everywhere. 
In terms of stability the behaviour is that of short beams in our denomination. The upper boundary condition, here 
clamping, plays a crucial role, and the instabilty may be divergence or flutter. If the length is much larger than Lg, 
so that £ ^ 1, gravity effects dominate in all the upper part. This results in a confinement of the instabilty in the 
lower part of the beam, as it appeared in all the mode shapes. The upper boundary condition plays a lesser role. 
It is remarkable that the characteristics of the instability of these long beams seem to be much simpler than that of 
short beams: flutter prevails in all cases, and the critical load, the frequency and the mode shape have very simple 
evolutions. 

For long beams flutter is reached when the load q is of order 1. In that range of loading the neutral point where 
tension vanishes, at x = g, is located near the limit between the two regions, x = 1. The critical load is then such 
that, in dimensionless variables 


P ~ (45) 

which shows the balance between the loading P and the two stiffnesses of the systems, El for the flexural rigidity 
and mg for the tension induced by gravity. 

It should be recalled here that damping was not considered in our calculations. Damping may significantly modify 
the critical value for flutter instability even when small 0, and has been found to have both a stabilizing and 
destabilising effects. In some cases, such as that of Beck’s column, it has been found that addition of damping 
destabilizes the system and changes the instability type from flutter to buckling. In gyroscopic systems such as fluid 
conveying-pipes or plates in axial flows, damping was found to destabilise neutral waves at zero frequency , which 
may also cause the instability to change from flutter to buckling in finite length systems. The generalisation to the 






11 


present case of a semi-infinite system tensionned by gravity needs to be addressed to see how damping could alter the 
buckcling/fiutter stability maps, and if it prevents flutter to prevail for all values of the follower force coefficient 77 . 


B. Relation to flow-induced instabilities 


The equivalence between the loading produced by the flow along a structure and a follower force exerted at its 
free end is discussed in [mil- Depending on the geometry of the problem, the load maybe either a fully follower 
load (fluid-conveying pipe) or partially a follower load (cylinder with axial flow). In the latter case, the shape of the 
downstream free end strongly influences the value of rj. 

For systems mainly tensionned by gravity, such as hanging fluid-conveying pipes or hanging ribbons under axial 
flow [71, [l2| the equations of motion can be put in the common dimensionless form 


dx* 


A 

dx 


(2 




(46) 


where the dimensionless variables have been defined using the length Lg, as we did in the present paper. Here v 
is a dimensionless flow velocity and /3 expresses the proportion of the fluid mass in the total mass. The boundary 
conditions at the lower end read 




(47) 


which expresses the case of a purely follower load. Eq. (HSl) is identical to (flUl) , except for the gyroscopic term which 
is dependent on /3. The stability diagram for the case of a long hanging pipe has been analysed by Q and for long 
hanging ribbons by [II- The critical velocity Vc was found to depend on the length in a manner very similar to 
that found in this paper. For very long pipes flutter arises at a critical velocity Vc that depends only on j3. Defining 
Qc = , and in the particular limit of /3 = 0, the results of Q should converge to that found in this paper for 77 = 1. 

This is not exactly the case when the results of Q are used, as there calculations have been done for small but non 
vanishing values of B. It is a well known feature that the case /? = 0 is only found as the limit of very small values 

( 10 - 3 ) oiBm- 

If the tensioning load is not gravity but a flow-induced friction on the slender structure, the form of the problem 
is changed in several ways [l^. First, the tensioning load increases with the flow velocity, which does not allow 
one to define an equivalent to the length Lg, as for gravity. The reference length that can then be defined scales 
the flow-induced forces that are proportionnal to volume (added stiffness forces) and those that are proportionnal 
to surface (friction). Second, the friction load is not constant in direction, as is gravity, but acts tangentially to 
the instantaneous position of the beam. Third, in the case of axial flow outside of a structure, the flow may also 
induce forces perpendicular to the instantaneous position of the beam, in the form of a transverse drag. Moreover, 
the boundary condition at the free end include other terms that depend on the local geometry of the flow. If all time 
dependent terms are neglected, to determine the static stability of the beam under flow, the equation defining the 
deflection reads [3 


as 2 d 




= 0 . 


(48) 


Upon defining x = and q = , the equation for the static behaviour of the beam, Eq. (II3, is recovered. If 

the end-condition sustains a force only in the beam axis, like in reference [l3l| . the results at 77 = 0 are recovered. For 
the analysis of the dynamic instability not only inertial terms need to be considered, as in the case of this paper. A 
gyroscopic term, as in Eq. ()46|) . and a damping term appear, further complicating the behaviour of the system. In 
the analysis of long towed cylinder Q edge flutter was also found, confirming the generality of the results presented 
here in a generic case. 
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